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We calculate the finite temperature and non-equilibrium electric current through systems de- 
scribed generically at low energy by a singlet and two spin doublets for TV and N ± 1 electrons 
respectively, coupled asymmetrically to two conducting leads, which allows for destructive interfer- 
ence in the conductance. The model is suitable for studying transport in a great variety of systems 
such us aromatic molecules, different geometries of quantum dots and rings with applied magnetic 
flux. As a consequence of the interplay between interference and Kondo effect, we find changes by 
several orders of magnitude in the values of the conductance and its temperature dependence as 
the doublet level splitting is changed by some external parameter. The differential conductance at 
finite bias is negative for some parameters. 

PACS numbers: 73.23.-b,73.22.-f, 75. 20. fir 



Transport properties of single molecules are being ex- 
tensively studied due to their potential use as active com- 
ponents of new electronic devices. Experimental results 
show an increased conductance at low temperatures due 
to the usual spin-1/2 or spin-1 Kondo effect and quantum 
phase transitions were induced changing externally con- 
trolled parameters [H-[H|- Measurements through single 
7T-conjugated molecules [J, stimulated further theo- 
retical work. In particular, the possibility of controlling 
molecular electronics using quantum interference effects 
in annulene molecules has been recently proposed 
The conductance depends on which sites of the molecule 
are coupled to the leads and on interference phenomena 
related to the symmetry of the system. For certain con- 
ditions it can be totally suppressed and restored again 
by symmetry-breaking perturbations [l2| . 



Interference phenomena are also well known in rings 
threaded by a magnetic flux. While sizable fluxes are at 
present impossible to apply to small annulene molecules 
due to their small area, Aharonov-Bohm oscillations 
were observed in systems involving two quantum dots 
(QDs).[H Systems of three [30], and more Q Ds 
have been assembled to study the effects of interdot hop- 
ping on the Kondo effect, and other physical properties 
driven by strong correlations. It has been predicted that 
the transmittance integrated over a finite energy window 
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u or the conductance through a ring of strongly cor- 
related one-dimensional systems 2(| 2lJ display dips as 
a function of the applied magnetic flux at fractional val- 
ues of the flux quantum, due to spin-charge separation. 
While the energy integration mimics a finite bias voltage 
V or temperature T applied to the system, the calcula- 
tions were actually performed at V = T = 0. Moreover, 
as in calculations of transport through molecules which 




FIG. 1. Scheme of the relevant matrix elements at low ener- 
gies. 



were included perturbatively, missing the Kondo regime, 
for which the conductance is usually the highest [23| . 
The Kondo effect was also missed in previous studies of 
interference effects at equilibrium (V = 0) for three dots 
at T > Tk [241 ] . where Tk is the Kondo temperature of 
the system, and in the spinless case (25|. Experimen- 
tally, the crossing of levels and the ensuing destructive 
interference has been induced applying magnetic field to 
a system with large, level-dependent g factors (26|. 

An analysis of the current through finite rings shows 
that interference phenomena take place when two lev- 
els become degenerate 
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included the effect of correlations [llj, |12|, [22j , the leads 



molecules these states correspond to two doublets with 
total wave vectors ±K , degenerate due to reflection sym- 
metry (in absence of an external flux) [0 , 12 , 21 1 . For the 
infinite-?/ Hubbard model, the interfering states differ in 
the spin quantum numbers [2lj [ . 

In this work we study the non-equilibrium transport in 
a generalized Anderson model, whose localized part con- 
sists of a singlet with even N electrons and two doublets 
with iV±l electrons (see Fig. [1]). The model describes the 
low-energy physics and interference phenomena of a wide 
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range of physical systems, such as those mentioned above. 
This calculation is a substantial improvement to previous 
results on conductance through interacting systems with 
interfering effects, 12, IH-21| since it introduces the rel- 
evant Kondo regime neglected in those studies as well as 
finite temperature and non-equilibrium effects at finite 
gate voltages. 

Due to its reliability, for this study we use the non- 
crossing approximation (NCA) [27, 28 1. For the case of 
one doublet, comparison of NCA with numerical renor- 
malization group (NRG) results [13] , shows that the NCA 
describes accurately the Kondo physics and in addition, 
allows us to reach finite bias voltages and temperatures. 
The leading behavior of the differential conductance with 
voltage is the correct one and G(T)/G(0) practically co- 
incides with the NRG result over several decades of tem- 
perature [29( . The results should improve with increasing 
degeneracy (small or zero level splitting in our case) and 
in fact experimental observations near quantum phase 
transitions were well reproduced by the NCA 0,|2^|. We 
calculate the effects of interference on non-equilibrium 
conductance in a regime of gate voltages for which the 
doublets are favored, leading to Kondo effect and high 
values of the conductance. In this approximation we 
asume that other higher-lying levels are negligible for the 
range of parameters considered. 

The effective Hamiltonian is 
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where the singlet |0) and the two doublets \ia) (i = 1, 2; 
a =t or I) denote the localized states, c\ kty create con- 
duction states in the left (v = L) or right (v = R) 
lead, and VJ describe the hopping elements between the 
leads and both doublets, assumed independent of k. By 
a gauge transformation, three of the four V£ can be 
made real and positive. If the system is a ring with 
the leads v connected at sites j u , the remaining effec- 
tive hopping can be chosen as Vr — |V^|e -8 ^, where 
6 = (ifl— Jl)(K2 — K\), and Ki is the wave vector of \ia) 
[2l|. The magnitude of the hoppings is given by |VJ| = 
|i„(icr|d]^ (T |0)|, where t v is the hopping between site j v 
and the lead v, and eft creates an electron (or a hole 
depending on the sign of the gate voltage) at site j [3lj |. 
If both doublets are related by a symmetry operation in 
the absence of an applied magnetic flux (for example re- 
flection symmetry), then \V£\ = \Vj?\. If, in addition, 
<j> = 7T and E 2 = Ei, the state \Lo) = (\la) + \2o))/y/2 
(\Ra) = (|lcr) — |2cr))/>/2) mixes only with the left (right) 
lead and therefore current cannot be transported between 
the leads. 

As a basis for our study, we start from this situation of 
perfect destructive interference ((f) = tt) and also assume 



for simplicity symmetric coupling to the leads, so that 
\V$\ — V'/v2 and allow for a finite splitting of the dou- 
blets 5 = E2 — Ei (introduced for example by an applied 
magnetic flux [21, 24]). For 6 = 0, Eq. (JJ takes the form 
of an SU(4) Anderson model (with on-site hybridization 
V), which was used to interpret transport experiments in 



carbon nanotubes 32j. However, the system represented 



is different and the conductance vanishes in our case. An 
important technical difference is that in our model, the 
hybridization matrices of the states \ia) with the left and 
right leads are not proportional for </> ^ 0, and as a con- 
sequence tricks used to relate the conductance at V = 
with the spectral density of states cannot be used and 
we have to calculate the conductance G = dl/dV by nu- 
merical differentiation of the current even for V — > 0, 
using a non equilibrium formalism [33| . Changing basis 
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(p}, it is seen that 6 acts as a symmetry breaking field 
on the SU(4) Anderson model, reducing the symmetry to 
SU(2). For d —> 00, the doublet with energy E2 can be 
neglected and the model reduces to the usual one-level 
SU(2) Anderson model. Therefore, the model interpo- 
lates between the one-level SU(4) and SU(2) Anderson 
models. 

At T = V — 0, the conductance G(T, V) is given in 
terms of the scattering phase shifts in a Fermi liquid de- 
scription [3^. In turn, these phase shifts can be related 
to the expectation values rii a = (\ia)(ia\) generalizing 
the Friedel sum rule t 35] to the SU(4) model with a sym- 
metry breaking field. For constant density of states of 
the leads we obtain 



G = G(0, 0) = — V sin 2 
h ^-^ 
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For the numerical calculations, we assume a constant 
density of states per spin of the leads p between — D 
and D and take T = 2A = 2irp(V ) 2 as the unit of en- 
ergy. Without loss of generality, we take the Fermi level 
ep = E s =0, and assume 8 = E 2 — E\ > 0. We also 
limit our present study to gate voltages that drive the 
system to the Kondo regime E s — E\ 3> A, for which 
the conductance is highest. In Fig. [5] we compare Eq. 
@ as a function of 5, with the result obtained by nu- 
merical differentiation of the current calculated at very 
low temperatures (specifically T = 0.05Tk{S), see be- 
low). The agreement is very good and gives confidence 
on the numerical procedure and on the consistency of 
the NCA results for the current and occupation numbers 
rii = 2rii a , also displayed in Fig. [2] We define the Kondo 
temperature Tk as the half width at half maximum of the 
peak nearest to the Fermi energy of the spectral density 
pi a (uj) of the lowest doublet. As it is apparent from the 
inset in Fig. [21 we find that to a high degree of accuracy 
Tk = fTft, where / is a factor of the order of 1 (0.606 
for the parameters of Fig. [5J), and TY- is given by the 
following expression 
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FIG. 2. Equilibrium conductance and T <C Tk (left scale) 
and occupation numbers (right scale) as a function of the level 
splitting. Squares were obtained by numerical differentiation 
of the current and circles correspond to Eq. ((2j). Inset: Kondo 
temperature as a function of 8. Squares (solid line) correspond 
to the NCA (analytical) result. Parameters are A = 0.5. 
D = 10, Ei = -4, E 2 =Ei + S. 



Tk = {(D + S)Dexp [tt£i/(2A)] + <5 2 /4} 1/2 - 5/2, (3) 

obtained minimizing the energy of the simple variational 
wave function 

|V) =a\s) +J2M^)(0hka\S), 

ika 

where \S) is the many-body singlet state with the filled 
Fermi sea of conduction electrons and the state |0) at the 
localized site. 

Eq. ([3]) interpolates between both one-level SU (N) lim- 
its, T% = L>exp [nE 1 /(NA)} (N = 4 for 5 = and N = 2 
for i5 — > +oo). Roughly, Tx{5) remains constant at the 
SU(4) value T K (0) (0.016 for the parameters of Fig. |2J) 
as long as 5 < Tk(0), and then decreases nearly expo- 
nentially (by almost two orders of magnitude for 5 = 1) 
before flattening at the one-level SU(2) value. As it is 
apparent in Fig. [21 Tr-(O) is also the characteristic en- 
ergy scale for the variation of the conductance with 5, 
for T = V = 0. For 5 one order of magnitude less than 
T K (0), G is very small, while for 5 > T K (0), G ap- 
proaches the ideal value for one SU(2) doublet, 2e 2 /h. 

In Fig. [31 we show the evolution with temperature 
of the conductance for several values of the level split- 
ting 5. Note that the ordinate axis is normalized by Go 
and the abscissa by T^-, the value of the temperature for 
which G(T) = G(0)/2. As expected Tg- - T K . Both G 
and T£ have a strong variation with 5. Near the SU(4) 
limit 5 — > 0, Go is much smaller, but Tk is much larger; 
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FIG. 3. Equilibrium conductance as a function of temper- 
ature for several level splittings. Dashed [dashed-dot] line 
corresponds to the limit of one-level SU(4) [SU(2)] Ander- 
son model. Other parameters as in Fig. [2] For 5 = 0.0035, 
Go = 0.08 and T* K = 0.0137, while for S = 0.06, G = 0.97 
and T* K = 0.0025. 



therefore G(T) persists nearly constant for a much wider 
range of temperatures. Instead, in the one-level SU(2) 
limit of large 5, G(T) starts near ideal values but decays 
much faster with temperature. Near the SU(4) limit, 
the temperature dependence of the conductance is quite 
similar to that calculated for the simpler geometry of 
carbon nanotubes with one particle in the quantum dot 
[3^ , although in the present case, the magnitude is much 
smaller. The bump observed in the curve for 5 = 0.06 is 
due to the contribution of the excited doublet at energy 
E2. 
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FIG. 4. Differential conductance as a function of bias voltage 
at different temperatures for 5 = 0.0112. Other parameters 
as in Fig. [2] 

In Fig. 21 we show the non-equilibrium differential con- 
ductance as a function of bias voltage V for a level split- 
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FIG. 5. Current as a function of level splitting 5 for eV = 
T% — 0.016 and different temperatures proportional to Tk(S). 
Other parameters as in Fig. [2] 



ting 8 that corresponds to an intermediate region between 
the SU(4) and one-level SU(2) limits (5 ~ Tr- ~ 0.01), 
and for several temperatures. One observes negative 
G(V) for voltages above the characteristic energy scale 
5 ~ Tk- Our interpretation of this result is the follow- 
ing. Since 6 ~ Tk, there is only a partial effect of the 
interference at small bias voltages, and the lowest lying 
doublet plays a dominant role. However as the energy of 
the bias voltage eV increases beyond the level splitting 6, 
both doublets contribute with nearly equal weight to the 
current, but in opposite ways due to destructive inter- 
ference. Therefore, the current starts to decrease as the 
effect of the excited doublet increases for larger applied 
bias voltages. 

In Fig. [SJ we display the current as a function of 5 for 
a fixed bias voltage at low temperatures. The bias volt- 
age was chosen to correspond to the characteristic energy 
Tr(0) of the SU(4) limit. Remarkably, a non monotonic 
behavior of the current is obtained due to two competing 
effects. One would expect a monotonous increase of the 
current with the level splitting as a consequence of the 
weakening of the destructive interference. However, the 
Kondo temperature Tk decreases strongly with 6, and for 
a fixed bias voltage a reduction of Tk implies a decrease 
in the current: at low V and T = 0, for one level with 
SU(N) symmetry, G(V) = dl/dV has a peak of width of 
the order of 2Tk/c Therefore the current at bias volt- 
ages which exceed a few Tk /e is roughly proportional to 
T K . 

We have studied the finite temperature and non equi- 
librium transport properties of an effective impurity 
Anderson model containing two doublets, which de- 
scribes the low-energy physics of quantum interference 
in nanoscopic systems. In the Kondo regime, dramatic 



changes in the values of the conductance and its temper- 
ature dependence take place as the doublet level split- 
ting is changed by some external parameter. For total 
destructive interference, the model interpolates between 
the SU(4) Anderson model when the splitting of the two 
doublets is S = 0, and the usual SU(2) model for large 
6. In the Kondo regime, while the characteristic temper- 
ature Tk increases significantly towards the SU(4) limit, 
both, the equilibrium conductance (gate voltage V = 0) 
and the total current at finite gate voltages vanish due 
to destructive interference. For finite V the total current 
peaks at finite i5 due to the interplay between interfer- 
ence and Kondo effects and when 5 ~ Tk, the differential 
conductance becomes negative due to partial destructive 
interference. 

In summary, by considering the interplay between two 
relevant effects, i.e. quantum interference and Kondo 
screening, we have shown the important consequences 
this can have on transport properties through a great 
variety of nanoscopic and molecular systems. 
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